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We investigate the general relativistic dynamics of Robertson-Walker models with a non-linear 
eqnation of state (EoS), focusing on the quadratic case P — Pq + ap -\- fip^ ■ This may be taken to 
represent the Taylor expansion of any arbitrary barotropic EoS, P{p)- With the right combination 
of Po, Q and (3, it serves as a simple phenomenological model for dark energy, or even unified 
dark matter. Indeed we show that this simple model for the EoS can prodnce a large variety of 
qualitatively different dynamical behaviors that we classify using dynamical systems theory. An 
almost universal feature is that accelerated expansion phases are mostly natural for these non-linear 
EoS’s. These are often asymptotically de Sitter thanks to the appearance of an effective cosmological 
constant. Other interesting possibilities that arise from the quadratic EoS are closed models that 
can oscillate with no singularity, models that bounce between infinite contraction/expansion and 
models which evolve from a phantom phase, asymptotically approaching a de Sitter phase instead 
of evolving to a “Big Rip”. In a second paper we investigate the effects of the quadratic EoS in 
inhomogeneous and anisotropic models, focusing in particular on singularities. 

PACS numbers: 98.80.Jk, 98.80.-k, 95.35.-|-d, 95.36.-l-x 


I. INTRODUCTION 


Model building in cosmology requires two main ingre¬ 
dients: a theory of gravity and a description of the mat¬ 
ter content of the universe. In general relativity (GR) 
the gravity sector of the theory is completely fixed, there 
are no free parameters. The matter sector is represented 
in the field equations by the energy-momentum tensor, 
and for a fluid the further specification of an equation of 
state (EoS) is required. Apart from scalar fields, typical 
cosmological fluids such as radiation or cold dark matter 
(CDM) are represented by a linear EoS, P = wp. 

The combination of cosmic microwave background ra¬ 
diation (CMBR) Q, 0. large scale structure (LSS) 0 
and supernova type la (SNIa) 0 observations provides 
support for a flat universe presently dominated by a com¬ 
ponent, dubbed in general “dark energy”, causing an ac¬ 
celerated expansion. The simplest form of dark energy 
is an ad hoc cosmological constant A term in the field 
equations, what Einstein called his “biggest blunder”. 
However, although the standard ACDM “concordance” 
model provides a rather robust framework for the inter¬ 
pretation of present observations (see e.g. 0S) , it re¬ 
quires a A term that is at odds by many order of magni¬ 
tudes with theoretical predictions 0 . This has prompted 
theorists to explore possible dark energy sources for the 
acceleration that go beyond the standard but unsatisfac¬ 
tory A. With various motivations, many authors have 
attempted to describe dark energy as quintessence, k- 
essence or a ghost field, i.e. with scalar fields with various 
properties. There have also been attempts to describe 
dark energy by a fluid with a specific non-linear EoS like 
the Chaplygin gas 0 , generalized Chaplygin gas 0 , van 
der Waals fluid 0 , wet dark fluid ^ and other spe¬ 


cific gas EoS’s [T3- Recently, various “phantom models” 
(ly = Pjp < —1) have also been considered min 
More simply, but also with a higher degree of generality, 
many authors have focused on phenomenological mod¬ 
els where dark energy is parameterized by assuming a 
w = P/p = wia). where a = a{t) is the expansion scale 
factor (see e.g. |l^ll^ ). 

Another possibility is to advocate a modified the¬ 
ory of gravity. At high energies, modification of grav¬ 
ity beyond general relativity could come from extra di- 
mensions, as req uired in string theory. In the brane 
world [ 13 , in, m m scenario the extra dimensions pro¬ 
duce a term quadratic in the energy density in the effec¬ 
tive 4-dimensional energy-momentum tensor. Under the 
reasonable assumption of neglecting 5-dimensional Weyl 
tensor contributions on the brane, this quadratic term 
has the very interesting effect of suppressing anisotropy 
at early enough times. In the case of a Bianchi I brane- 
world cosmology containing a scalar field with a large 
kinetic term the initial expansion is quasi-isotropic |2l| . 
Under the same assumptions, Bianchi I and Bianchi 
V brane-world cosmological models containing standard 
cosmological fluids with linear EoS also behave in a sim¬ 
ilar fashion^ |^, and the same remains true for more 
general homogeneous models |2^l^ and even some in¬ 
homogeneous exact solutions |^. Finally, within the 
limitations of a perturbative treatment, the quadratic- 
term-dominated isotropic brane-world models have been 


^ This only requires P/p = ro > 0, as opposed to in > 1 in the GR 
case. In the case of ekpyrotic/cyclic and pre-big bang models the 
init ial expansion is only isotropic if ro > 1 as in the case of GR 

0 . 
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shown to be local past attractors in the larger phase space 
of inhomogeneous and anisotropic models • More 

precisely, again assuming that the 5-d Weyl tensor con¬ 
tribution to the brane can be neglected, perturbations of 
the isotropic models decay in the past. Thus in the brane 
scenario the observed high isotropy of the universe is the 
natural outcome of generic initial conditions, unlike in 
GR where in general cosmological models with a stan¬ 
dard energy momentum tensor are highly anisotropic in 
the past (see e.g. H^). 

Recently it has been shown that loop quantum gr avity 
corrections result in a modihed Friedmann equation pUTf , 
with the modihcation appearing as a negative term which 
is quadratic in the energy density. Further motivation 
for considering a quadratic equation of state comes from 
recent studies of fc-essence fields as unified dark matter 
(UDM) models^ 113 • The general fc-essence field can 
be described by a fluid with a closed-form barotropic 
equation of state. The UDM fluid discussed in [sj has a 
non-linear EoS of the form P oc at late times. More 
recently, it has been shown that any purely kinetic k- 
essence field can be interpreted as an isentropic perfect 
fluid with an EoS of the form P = P{p) [13. Also, low 
energy dynamics of the Higgs phase for gravity have been 
shown to be equivalent to the irrotational flow of a per¬ 
fect fluid with equation of state P = Ill- 

Given the isotropizing effect that the quadratic energy 
density term has at early times in the brane scenario this 
then prompts the question: can a term quadratic in the 
energy density have the same effect in general relativity. 
This question is non-trivial as the form of the equations 
in the two cases is quite different. On the brane, for 
a given EoS the effective 4-dimensional Friedmann and 
Raychaudhuri equations are modihed, while the continu¬ 
ity equation is identical to that of GR. With the intro¬ 
duction of a quadratic EoS in GR, the Friedman equation 
remains the same, while the continuity and Raychaudhuri 
equations are modihed^. 

Taking into account this question (to be explored in 
detail in Paper II 0), the diverse motivations for a 
quadratic energy density term mentioned above and with 
the dark energy problem in mind, in this paper we explore 
the GR dynamics of homogeneous isotropic Robertson- 
Walker models with a quadratic EoS, P = Po + ap + Pp"^. 
This is the simplest model we can consider without mak¬ 
ing any more specihc assumptions on the EoS (H. It 
represents the hrst terms of the Taylor expansion of any 
EoS function P = P{p) about p = 0. It can also be 
taken to represent (after re-grouping of terms) the Tay¬ 
lor expansion about the present energy density po, see 
[13 . In this sense therefore the out-coming dynamics is 


^ These attempt to provide a unified model for both the dark mat¬ 
ter and the dark energy components necessary to make sense of 
observations. 

® With Respect to the case of the same EoS with vanishing 
quadratic term. 


very general. Indeed it turns out that this simple model 
can produce a large variety of qualitatively different dy¬ 
namical behaviors that we classify using dynamical sys¬ 
tems theory [13 [13 ■ An outcome of our analysis is that 
accelerated expansion phases are mostly natural for non¬ 
linear EoS’s. These are in general asymptotically de Sit¬ 
ter thanks to the appearance of an effective cosmological 
constant. This suggests that an EoS with the right com¬ 
bination of Pq) cr and /3 may provide a good and simple 
phenomenological model for UDM, or at least for a dark 
energy component. Other interesting possibilities that 
arise from the quadratic EoS are closed models that can 
oscillate with no singularity, models that bounce between 
infinite contraction/expansion and models which evolve 
from a phantom phase, asymptotically approaching a de 
Sitter phase instead of evolving to a “big rip” or other 
pathological future states 0, [33, [13 ■ 

As mentioned before, the question of the dynamical 
effects the quadratic energy density term has on the 
anisotropy in GR is explored in Paper II [13. There 
we analyze Bianchi I and V models with the EoS P — 
ap + Pp^ , as well as perturbations of the isotropic past 
attractor of those models that are singular in the past. 
We anticipate that Bianchi I and V non-phantom mod¬ 
els with /3 > 0 have an isotropic singularity, i.e. they 
are asymptotic in the past to a certain isotropic model, 
and that perturbations of this model decay in the past. 
Phantom anisotropic models with /3 > 0 are necessar¬ 
ily asymptotically de Sitter in the future, but the shear 
anisotropy dominates in the past. For /3 < 0 all mod¬ 
els are anisotropic in the past, while their specihc future 
evolution depends on the value of a. 

The paper is organized as follows. In section im we 
outline the setup and the three main cases we will inves¬ 
tigate. In section Unl we study the dynamics of isotropic 
cosmological models in the high energy limit (neglect¬ 
ing the Pq term). We hnd the critical points, their sta¬ 
bility nature and the occurrence of bifurcations of the 
dynamical system. In section nvi we consider the low 
energy limit (neglecting the p^ term). The full system 
is then analyzed in section El showing the qualitatively 
different behavior with respect to the previous cases. We 
then hnish with some concluding remarks and an outline 
of work in progress in section ED Units are such that 
SttG/c^ = 1. 


II. COSMOLOGY WITH A QUADRATIC EOS 

A. Dynamics with non-linear EoS 

The evolution of Robertson-Walker isotropic models 
with no cosmological constant A term is given in GR by 
the following non-linear planar autonomous dynamical 
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system: 

p = -ZH{p + P), (1) 

H = -H^-\{p + iP), (2) 

where P[ is the Hubble expansion function, related to 
the scale factor a hy H = d/a. In order to close this 
system of equations, an EoS must be specified, relating 
the isotropic pressure P and the energy density p. When 
an EoS P = P{p) is given, the above system admits a 
first integral, the Friedman equation 

H^ = \p-^, (3) 

3 a‘^ 

where K is the curvature, K = 0,±1 as usual for flat, 
closed and open models. 

Here we are interested in exploring the general dynami¬ 
cal features of a non-linear EoS P = P{p). Before consid¬ 
ering the specific case of a quadratic EoS, we note some 
important general points. 

First, it is immediately clear from Eq. o that an ef¬ 
fective cosmological constant is achieved whenever there 
is an energy density value pA such that P{pa) = —PA- 
More specifically: 

Remark 1. If for a given EoS function P = P{p) there 
exists a pA such that P{pa) = — Pa, then pA has the 
dynamical role of an effective cosmological constant. 

Remark 2. A given EoS P{p) may admit more than 
one point pA- If these points exist, they are fixed points 
of Eq. (Q. 

Remark 3. From Eq. since H -|- = a/o, an 

accelerated phase is achieved whenever P{p) < —p/3. 

Remark 4. Remark 3 is only valid in GR, and a 
different condition will be valid in other theories of 
gravity. Remarks 1 and 2, however, are only based 
on conservation of energy, Eq. O- The latter is also 
valid (locally) in inhomogeneous models, provided that 
the time derivative is taken to represent the derivative 
along the fluid flow lines (e.g. see ^ direct 

consequence of = 0. Thus Remarks 1 and 2 are 
valid in any gravity theory where = 0, as well as 

(locally) in inhomogeneous models. 

Second, assuming expansion, H > 0, we may rewrite 
Eq. (Q as: 

g = -3[p + P(p)], (4) 

where t = In a. Eq. 0 is a 1-dimensional dynamical 
system with fixed point(s) pa(s), if they exist. If 
p + P(p ) < 0 the fluid violates the null energy condi¬ 
tion |4l|, [431 and E g. (fTll im plie s what has been dubbed 
phantom behavior (cf. |43), i.e. the fluid behaves 


counter intuitively in that the energy density increases 
(decreases) in the future for an expanding (contracting) 
universe. Then: 

Remark 5. Any point pA is an attractor (repeller) of 
the evolution during expansion (the autonomous system 
0) if p -I- P{p) < 0 (> 0) for p < PA and p -I- P{p) > 0 
(< 0) for p > PA- 

Remark 6. Any point pA is a shunt^ of the autonomous 
system Eq. 0 if either p + P{p) < 0 on both sides of pA, 
or p -I- P(p) > 0 on both sides of p > pA. In this case the 
fluid is respectively phantom or standard on both sides. 

Let’s now consider the specific case of a general 
quadratic EoS of the form: 

P = Po + ap + f3p^. (5) 

The parameter (3 sets the characteristic energy scale pc 
of the quadratic term as well as it’s sign e 

/?=-. ( 6 ) 

Pc 

Remark 7. Eq. © represents the Taylor expansion, 
up to 0(3), of any barotropic EoS function P = P{p) 
about p = 0. It also represents, after re-grouping of 
terms, the Taylor exp ansion about the present energy 
density value po In this sense, the dynamical 

system (tna with © is general, i.e. it represents 
the late evolution, in GR, of any cosmological model 
with non-linear barotropic EoS approximated by Eq. ©. 

The usual scenario for a cosmological fluid is a stan¬ 
dard linear EoS (Pq = /3 = 0), in which case a = ui is 
usually restricted to the range — 1 < a < 1. For the 
sake of generality, we will consider values of a outside 
this range, considering dynamics only restricted by the 
request that p > 0. The first term in Eq © is a con¬ 
stant pressure term which in general becomes important 
in what we call the low energy regime. The second term 
is the standard linear term usually considered, with 


If it is positive, a has an interpretation in terms of the 
speed of sound of the fluid in the limit p ^ 0, a = c^. 
The third term is quadratic in the energy density and will 
be important in what we call the high energy regime. 

In the following, we first split the analysis of the dy¬ 
namical system Eqs. (D m 0 in two parts, the high 
energy regime where we neglect Pq the low energy 


^ This is a fixed point which is an attr act or for one side of the 
phase line and a repeller for the other I 33 
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regime where we set /3 = 0, then we consider the full sys¬ 
tem with EoS ©■ Using only the energy conservation 
Eq. Q we list the various sub-cases, also briefly antic¬ 
ipating the main dynamical features coming out of the 
analysis in Sections HTTl IIVI and Ivl 


B. Quadratic EoS for the high energy regime 

In the high energy regime we consider the restricted 
equation of state: 


Phe = ap-\ -. (8) 

Pc 

The energy conservation Eq. o can be integrated in gen¬ 
eral to give: 


PHE{a) = 
A = 


A(a -h l)pc 
a3(“+i) - eA’ 

3(a+l) 

Poao 

{a + l)pc + epo ’ 


(9) 

( 10 ) 


where Po, Uo represent the energy density and scale factor 
at an arbitrary time to- This is valid for all values of e, 
Pc and a, except for a ^ —1. In the case a = —1 the 
evolution of the energy density is: 


PHE{a) 



Po Pc \ J 


( 11 ) 


The EoS with this particular choice of parameters has 
already been considered as a possible dark energy model 
0|. We will concentrate on the broader class of 
models where a ^ —1. 

In Section Hm we will give a dynamical system analysis 
of the high energy regime, but it is first useful to gain 
some insight directly from Eq. 0. 

We start by defining 


PA :=-e(l + a)pc, (12) 

noticing that this is an effective positive cosmological 
constant point only if e(l -I- a) <0. It is then convenient 
to rewrite Eq. m in three different ways, defining 
a* = |A|^/3(a-i-i) ^ each representing two different sub¬ 
cases. 


A: e(l -I- a) > 0, PA < 0, 


P = 


|1 + a\pc 
\3(l+a) 

:) 


(13) 


Al: e > 0, (1 -I- a) > 0. In this case a* < a < oo, with 
oo > p > 0. Further restrictions on the actual range of 
values that a and p can take may come from the geome¬ 
try. For a subset of appropriate initial conditions closed 
(positively curved) models may expand to a maximum a 


(minimum p) and re-collapse, and for a < —1/3 not all 
closed models have a past singularity at a = a*, having 
instead a bounce at a minimum a (maximum p). 

A2: e < 0, (1 -I- a) < 0. In this case 0 < a < a*, with 
0 < p < oo, and the fluid exhibits phantom behavior. 
All models have a future singularity at a = a*, but in 
general closed models contract from a past singularity, 
bounce at a minimum a and p, then re-expand to the 
future singularity (we will refer to this as a phantom 
bounce). 


B: PA > 0, p > PA, 


P = 


PA 


1 - 



3(l+a) ■ 


(14) 


Bl: e > 0, (1 -I- a) < 0, A > 0. In this case a* < a < oo, 
with oo > p > PA. As in case Al, further restrictions 
on the actual range of values that a and p can take 
may come from the geometry. For a subset of initial 
conditions closed models may expand to a maximum a 
(minimum p) and re-collapse, while for another subset 
closed models don’t have a past singularity at a = a*, 
having instead a bounce at a minimum a (maximum p). 


B2: e < 0, (1 -I- a) > 0, A < 0. In this case 0 < a < a*, 
with PA < p < oo. As in the case A2, the fluid has a 
phantom behavior. All models have a future singularity 
at a = a*, with closed models contracting from a past 
singularity to a minimum a and p before re-expanding. 


C: PA > 0, p < PA, 


P = 


PA 


1-b 



3(l+a) ■ 


(15) 


Cl: e > 0, (1 -I- a) < 0, A < 0. In this case 0 < a < oo, 
with 0 < p < p\. The fluid behaves in a phantom 
manner but avoids the future singularity and instead 
evolves to a constant energy density pA- Closed models, 
however, typically bounce with a minimum p at a finite a. 


C2: e < 0, (1 -I- a) > 0, A > 0. In this case 0 < o < oo, 
with PA > p > 0. Again, closed models may evolve 
within restricted ranges of a and p, even oscillating, for 
a > —1/3, between maxima and minima of a and p. 


C. Low energy regime: affine EoS 

In the low energy regime we consider the affine equa¬ 
tion of state: 


Ple = Po + ctp- (16) 

This particular EoS has been investigated as a possible 
dark energy model [ 11113 , however, only spatially flat 









Friedmann models where considered. The scale factor 
dependence of the energy density is: 


5 


PLE{a) = — 

B = 


Pn 


(a + 1) 

Po 


(l + a) 


3(l+c<) 


(17) 

(18) 


This is valid for all values of Po and a except a ^ —1. In 
the case a = — 1, the evolution of the energy density is: 


El: Po > 0, (1 + a) < 0, B > 0. In this case 0 < a < oo, 
with Pa < P < oo. As in the case D2, the fluid behaves 
in a phantom manner. The flat and open models are 
asymptotically de Sitter in the past, when their energy 
density approaches a finite value (p ^ pA as a ^ 0), 
and when pA becomes negligible in Eq. (PI they evolve 
as standard linear phantom models, reaching a future 
singularity in a finite time (p —> oo as a —> oo). The 
closed models contract to a minimum a (minimum p), 
bounce and re-expand. 


PLE{a) = Po-3Poln^^^, (19) 

As in the high energy case, we will concentrate on the 
broader class of models where a ^ —1. 

In Section ITvI we present the dynamical system analysis 
of the low energy regime, but first let us gain some insight 
from Eq. 1I3. As with the high energy case, in many 
cases the fluid violates the null energy condition (p + P < 
0) and exhibit phantom behavior. Defining 

Pa '■=—Po/(l + o:), (20) 

we see that a positive effective cosmological constant 
point exists, pA > 0, only if Po/(l + a) < 0. Eq. itTTI) 
can be rewritten in three different ways, defining 
a* = each representing two different sub¬ 

cases. 


D: Po/(l -I- a) > 0, PA < 0, 


Po 

(a + 1) 



( 21 ) 


Dl: Po >0, (1 -I- a) > 0. In this case 0 < a < oo, 
with oo > p > — |pa|- The geometry places further 
restrictions on the values that a and p can take. The 
subset of open models (negative curvature) are all 
non-physical as they evolve to the p < 0 region of the 
phase space. The spatially flat models expand to a 
maximum a (when p = 0) and recollapse. The closed 
(positively curved) models expand to a maximum a 
(minimum p) and recollapse, and for —I<a<—l/3a 
subset of closed models oscillate between a maximum 
and minimum a (minimum and maximum p). 


D2: Po < 0, (1 -I- a) < 0. In this case 0 < a < oo, 
with — |pa| < p < oo. In this case the fluid exhibits 
phantom behavior. The subset of open models are all 
non-physical as they evolve from the p < 0 region of the 
phase space. The spatially flat models contract, bounce 
at a minimum a when p = 0 and re-expand in the future. 
The closed models contract, bounce at a minimum a 
and p, then re-expand in the future. 


E: PA > 0, p > PA, 


E2: Po < 0, (1 -I- a) > 0, B > 0. In this case 0 < a < oo, 
with oo > p > PA. All flat and open models expand from 
a singularity and asymptotically evolve to a de Sitter 
model, with p = pA- The closed models evolve from a 
contracting de Sitter model to minimum a (maximum 
p), bounce and then evolve to an expanding de Sitter 
model. 


F: PA > 0, p < PA, 


P = PA- 



(23) 


FI: Po > 0, (1 -I- a) < 0, i? < 0. In this case 0 < a < oo, 
with PA > p > — oo. The subset of open models are 
all non-physical as they evolve to the p < 0 region 
of the phase space. The flat models evolve from an 
expanding de Sitter phase to a contracting de Sitter 
phase. The closed models oscillate between a maximum 
and minimum a (minimum and maximum p). 


F2: Po < 0, (1 -I- a) > 0, P < 0. In this case 0 < a < oo, 
with — oo < p < PA. The fluid exhibits phantom behav¬ 
ior. The open models are all non-physical as they evolve 
from the p < 0 region of the phase space. The flat and 
closed models evolve from a contracting de Sitter phase, 
bounce at minimum a and p, then re-expand, asymptot¬ 
ically approaching a expanding de Sitter phase. 


D. The full quadratic EoS 


In Section El we present the dynamical system analy¬ 
sis of the full quadratic EoS models given by Eq. (0 , but 
again we first study the form of p(a) implied by conser¬ 
vation of energy, Eq. O- As with the previous cases the 
fluid can violate the null energy condition (p+P < 0) and 
therefore may exhibit phantom behavior. The system 
may admit two (possibly negative) effective cosmological 
constant points: 


PA,i ■■= ^ -(a-f l)-f VA , 

(24) 

PA,2 ■■= ^ -(a-f 1)- VA , 

(25) 

A := (a -f 1)2 - Af3Po 

(26) 


P = PA + 


-3(l+a) 


( 22 ) 


if 
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is non negative. Clearly, the existence of the effective 
cosmological points depends on the values of the param¬ 
eters in the EoS. This in turn affects the functional form 
of p{a). In order to find p{a) the following integral must 
be evaluated: 


— 3 In 



dp 


Po -\- (pt l)p + 


(27) 


This is done separately for the cases when no effective 
cosmological points exist (A < 0), when one cosmo¬ 
logical point exist, pA,i = Pa, 2 = Pa 7^ 0 (A = 0) and 
when two cosmological points exist, ^ PA ,2 7 ^ 0 
(A >0). We now consider these three separate sub-cases. 


G: (1-b a)2 < 4/3Po, A < 0, 

r- v1^tan(|v^ln(^)) (a + 1 ) 

" ” + 2/3 ' 

r = 2Ppo + {a + \). (28) 

Gl: /? ^ 0, F^o d■ In this case ui u ^ 112 (ivhere 
tti < 02 ), with 00 > p > — 00 . The fluid behaves in a 
standard manner and all models have a past singularity 
at a = oi. All open models are non-physical as they 
evolve to the p < 0 region of the phase space. The 
flat models expand to a maximum a (p = 0) and then 
re-collapse. The closed models can behave in a similar 
manner to flat models except they reach a minimum 
p before re-collapsing. Some closed models oscillate 
between maxima and minima a and p. 


G2: /? < 0, Po < 0. In this case oi < a < 02 (where 
ai < 02 ), with — 00 < p < 00 . The fluid behaves in 
a phantom manner. All open models are non-physical 
as they evolve from the p < 0 region of the phase 
space. The flat and closed models represent phantom 
bounce models, that is they evolve from a singularity at 
a = ai (p = 00 ), contract to a minimum a (minimum 
p) and then re-expand to the future singularity at a = 02 . 


manner. If p\ > 0, the open and flat models evolve from 
a past singularity (a = oi) and evolve to a expanding 
de Sitter phase. For a subset of initial conditions closed 
models may expand to a maximum a (minimum p) 
and re-collapse, while for another subset closed models 
avoid a past singularity, instead having a bounce at a 
minimum a (maximum p). If p\ < 0, the open models 
are non-physical, while flat and closed models represent 
recollapse models. 


H3: P < 0, Po < 0, p < PA. In this case oi < a < 00 
with — 00 < p < Pa- The fluid behaves in a phantom 
manner. The open models are all non-physical as they 
evolve from the p < 0 region of the phase space. The flat 
and closed models evolve from a contracting de Sitter 
phase, bounce at minimum a and p, then re-expand, 
asymptotically approaching an expanding de Sitter 
phase. 


H4: p < 0, Po < 0, p > PA. In this case 0 < a < oi 
with Pa < p < 00 and the fluid behaves in a phantom 
manner. All models have a future singularity at a = oi. 
If PA > 0, closed models contract from a past singularity 
to a minimum a and p before re-expanding (phantom 
bounce), while flat and open models are asymptotic to 
generalized de Sitter models in the past. If pA < 0, open 
models are non-physical, while flat and closed models 
contract from a past singularity to a minimum a and p 
before re-expanding. 


I: {l + af > ApPo, A > 0, 


P 


PA ,2 



— PA,lC 


-c 


Po — PA,2 
Po — PA,1 


(30) 

(31) 


Note that /3 > 0 (< 0) implies pA ,2 < PA,i (pA,i < PA. 2 ), 
and G < 0 implies pA ,2 < Po < PA.i for P > Q 
(PA.i < Po < PA ,2 for P < 0). 


H: (1 + af = ApPo, A = 0, 


P = 



(29) 


HI: /3 > 0, Po > 0, p < PA- In this case 0 < o < oi with 
Pa > p > — 00 . The fluid behaves in a standard manner. 
The subset of open models are all non-physical as they 
evolve to the p < 0 region of the phase space. The flat 
models evolve from an expanding de Sitter phase to a 
contracting de Sitter phase. The closed models oscillate 
between maxima and minima a and p. 


H2: P > 0, Po > 0, p > PA- In this case oi < a < 00 
with 00 > p > PA and the fluid behaves in a standard 


II: P > 0, Po > 0, p < Pa, 2 ) hence we consider pA ,2 > 0. 
In this case 0 < a < ai with pA ,2 > p > —00 and the 
fluid behaves in a standard manner. The open models 
are all non-physical as they evolve to the p < 0 region 
of the phase space. The flat models evolve from an 
expanding de Sitter phase to a contracting de Sitter 
phase. The closed model region contains a generalized 
Einstein static fixed point and models which oscillate 
indefinitely (between minima and maxima a and p). 

12: P > 0, Po > 0, PA ,2 < p < PA.I- In this case 
0 < a < 00 with PA,2 < P < PA,i and the fluid behaves 
in a phantom manner. The open models evolve from 
one expanding de Sitter phase (p = pa,2) to more 
rapid (greater p and H) de Sitter phase (p = PA,i)j 
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however the spatial curvature is negative in the past 
and asymptotically approaches zero in the future. The 
flat models behave in a similar manner except that the 
curvature remains zero. The closed models undergo a 
phantom bounce with asymptotic de Sitter behavior, 
that is they evolve from a contracting de Sitter phase, 
reach a minimum o, minimum p and then evolve to a 
expanding de Sitter phase. 

13: /3 > 0, Po > 0, p > pA,i. In this case ai < o < oo 
with oo > p > pa,i and the fluid behaves in a standard 
manner. All flat and open models expand from a 
singularity at a = oi and asymptotically evolve to a 
expanding de Sitter phase {p = pa,i)- A subset of closed 
models evolve from a contracting de Sitter phase to 
minimum a (maximum p), bounce and then evolve to 
an expanding de Sitter phase. Another subset of closed 
models expand from a singularity at a = oi, reach a 
maximum a and minimum p, only to re-collapse. 

14: (3 < 0, Po < 0, p < pA,i- In this case ai < a < oo, 
with — oo < p < pa,i and the fluid behaves in a phantom 
manner. The open models are all non-physical as they 
evolve from the p < 0 region of the phase space. The flat 
and closed models evolve from a contracting de Sitter 
phase, bounce at minimum a and p, then re-expand, 
asymptotically approaching a expanding de Sitter phase. 


• “Big Bang”/“Big Crunch” : For a ^ 0, p —> oo. 

If the singularity occurs in the past (future) we refer to 
it as a “Big Bang” (“Big Crunch”). In order to differenti¬ 
ate between various types of singularities, we will use the 
following classification system for future singularities 
(cf. also 

• Type I (“Big Rip”) : For t ^ t*, a —> oo, p —> oo 
and |P| ^ oo. 

• Type II (“sudden”) : For t ^ t*, a ^ a*, p ^ p* 
or 0 and |P| ^ oo. 

• Type III : For t ^ p, a ^ a*, p —> oo and |P| ^ 
oo. 

• Type IV : For t —>■ ti,, a —> p —>■ p* or 0, |P| ^ 
|P*| or 0 and derivatives of H diverge. 

Here ti,, a*, p* and |P*| are constants with a* ^ 0. The 
main difference in our case is that the various types of sin¬ 
gularities may occur in the past or the future. The future 
singularity described in case A2 falls into the category of 
Type III, however, the past singularity mentioned in case 
AI is also a Type III singularity. In the case of the full 
quadratic EoS, all singularities which occur for a finite 
scale factor (a = oi) are of Type III. 


15: /3 < 0, Po < 0, pA,i < p < PA,2 (where pA,i < PA,2)- 
In this case 0 < a < oo with pA,2 > P > PA,i and the 
fluid behaves in a standard manner. The open models 
evolve from a expanding de Sitter phase (p = pa, 2) to 
less rapid (lower p and H) de Sitter phase (p = pA,i) 
with the spatial curvature being negative in the past 
and zero asymptotically in the future. The flat models 
behave in a similar manner, except that the curvature 
remains zero throughout the evolution. The closed 
models can undergo a phantom bounce with asymptotic 
de Sitter behavior in the future and past, a subset of 
these models enter a loitering phase both before and 
after the bounce. There are a subset of closed models 
which oscillate indefinitely. 

16: /3 < 0, Po < 0, p > PA,2- In this case 0 < a < oi 
with Pa, 2 < p < 00 and the fluid behaves in a phantom 
manner. All models have a future singularity at a = ai, 
with closed models contracting from a past singularity 
to a minimum a and p before re-expanding (phantom 
bounce). 


E. The Singularities 

In general, singularities may behave in qualitatively 
different ways. The singularities present for the non 
linear EoS are quite different from the standard “Big 
Bang”/“Big Crunch” singularity. The standard singu¬ 
larities are such that: 


III. HIGH ENERGY REGIME DYNAMIGS 

A. The dimensionless dynamical system 

It is convenient to describe the dynamics in terms of 
dimensionless variables. In the high energy regime these 
are: 


= -A : y = > V = '/W\t ■ 

\Pc\ VH 

The system of equations 0-0 then changes into: 
x' = —3y ((a + l)a: -I- ea;^) , 
y = -y'^ - ^ ((3a + l)x + 3ex'^) , 


(32) 


(33) 


and the Friedman equation © gives 


X K 

3 IpcIo^’ 


(34) 


The discrete parameter e denotes the sign of the 
quadratic term, e G { —1,1}. The primes denote differen¬ 
tiation with respect to 77 , the normalized time variable. 
The variable x is the normalized energy density and y the 
normalized Hubble function. We will only consider the 
region of the phase space for which the energy density re¬ 
mains positive (a; > 0). The system of equations above is 
of the form u( = fi{uj). Since this system is autonomous, 
trajectories in phase space connect the fixed/equilibrium 




points of the system (uj,o), which satisfy the system of 
equations fi{uj^o) = 0. The fixed points of the high en¬ 
ergy system and their existence conditions (the condi¬ 
tions for which x > 0 and x,y gM.) are given in Tabled 


TABLE I: Location and existence conditions (x > 0 and x,y £ 
R) of the fixed points of the high energy regime system. 


Name 

X 

y 

Existence 

M 

0 

0 

— CX3 < a < CX3 

E 

e(3a + l) 

3 

0 

e(3a -I- 1) < 0 

dS+ 

-h 

1 


e(a -I- 1) < 0 

dS- 

-h 

1 


e{a -1- 1) < 0 


The first fixed point (M) represents an empty flat 
(Minkowski) model. The parabola = x/3 is the 
union of trajectories representing flat models, K = 0 
in Eq. m (see Figs. 1 and 3-7). The trajectories be¬ 
low the parabola represent open models {K = — 1), while 
trajectories above the parabola represent closed models 
{K = -1-1). The second fixed point (E) represents a gen¬ 
eralized static Einstein universe. This requires some form 
of inflationary matter and therefore may only exists when 
a < —1/3 if e = -1-1 and when a > —1/3 if e = —1. 
The last two points represent expanding and contracting 
spatially flat de Sitter models (dS±). These points exist 
when the fluid permits an effective cosmological constant 
point, XA '■= Pa! Pc = —e(<a + l); in addition xa > 0 must 
be true for the fixed points to be in the physical region 
of the phase space. There are further fixed points at in¬ 
finity, these can be found by studying the corresponding 
compactified phase space. The first additional fixed point 
is at X = y = oo and represents a singularity with infinite 
expansion and infinite energy density. The second point 
isatx = oo,y = —oo and represents a singularity with 
infinite contraction and infinite energy density. 


B. Generalities of stability analysis 


The stability nature of the fixed points can be found 
by carrying out a linear stability analysis. In brief (see 
e.g. for details), this involves analyzing the behavior 
of linear perturbations Uj = Uj^o + Vj around the fixed 
points, which obey the equations = Mu^ . The matrix 
M is the Jacobian matrix of the dynamical system and 
is of the form: 


My 


dft 


duj 


Uk—Uk,o 


(35) 


The eigenvalues of the Jacobian matrix evaluated at 
the fixed points tell us the linear stability character of the 


fixed points. The fixed point is said to be hyperbolic if the 
real part of the eigenvalues is non-zero (K(Ai) ^ 0). If all 
the real parts of the eigenvalues are positive (R(Ai) > 0) 
the point is said to be a repeller. Any small deviations 
from this point will cause the system to move away from 
this state. If all the real parts are negative (K(Ai) < 0), 
the point is said to be an attractor. This is because if the 
system is perturbed away from this state, it will rapidly 
return to the equilibrium state. If some of the values 
are positive, while others are negative then the point is 
said to be a saddle point. If the eigenvalues of the fixed 
point are purely imaginary then the point is a center. 
If the center nature of the fixed point is confirmed by 
some non-linear analysis, then the trajectories will form 
a set of concentric closed loops around the point. If the 
eigenvalues do not fall into these categories, we will resort 
to numerical methods to determine their stability. 

The eigenvalues for the fixed points of the system 
(Eq.’s (123)) are given in Table QH and the linear stability 
character is given in Table cni 


TABLE II: Eigenvalues for the fixed points of the high energy 
regime system. 


Name 

Ai 

A2 

M 

0 

0 

E 


77 (^+ 1 ) 

dS+ 

(a - 1 - 1 )-^—3e(a -|- 1 ) 

— |7“3e(a -I- 1) 

dS- 

— (q; “h 1 )-^— 3e(Q; + 1 ) 

|7“3e(a - 1 - 1 ) 


TABLE III: The linear stability of the fixed points for the 
high energy regime system. 

Name 

e = +l 

e = -1 

M 

undefined 

undefined 

E 

Saddle (a 7 —1/3) 

Center (a 7 —1/3) 

dS+ 

Attractor (a < —1) 

Saddle (a > — 1) 

dS- 

Repeller (a < — 1) 

Saddle (a > — 1) 


C. The £ = -1-1 case 

We first consider the system when we have a posi¬ 
tive quadratic energy density term (e = -|-1) in the high 
energy regime EoS. We will concentrate on the region 
around the origin as this is where the finite energy den¬ 
sity fixed points are all located. The plots have been cre¬ 
ated using the symbolic mathematics application Maple 
9.5. The individual plots are made up by three layers, 
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FIG. 1: The phase space for the high energy regime system 
with e = +1 and a < — 1. The npper (lower) region corre¬ 
sponds to the case B1 (Cl). 


the first is a directional (represented by grey arrows) field 
plot of the state space. The second layer represents the 
separatrices and fixed points of the state space . A sep- 
aratrix (black lines) is a union of trajectories that marks 
a boundary between subsets of trajectories with differ¬ 
ent properties and can not be crossed. The fixed points 
are represented by black dots. The final layer represents 
some example trajectories (grey lines) which have been 
calculated by numerically integrating the system of equa¬ 
tions for a set of initial conditions. The character of the 
fixed point M is undefined and so is determined numer¬ 
ically. The fixed point representing the generalized Ein¬ 
stein static solution is a saddle point. The fixed points 
representing the generalized expanding (contracting) de 
Sitter points are attractor (repeller) points. The trajec¬ 
tories or fixed points in the y > 0 (y < 0 ) region represent 
expanding (contracting) models. We will mainly discuss 
the right hand side of the state space (expanding mod¬ 
els) as in general the corresponding trajectory on the left 
hand side is identical under time reversal. 


1. The a < —1 sub-case 

The phase space of the system is considered when 
a < — 1 and is shown in Fig. 1. The lowest horizontal 
line {x = 0) is the separatrix for open models {K = —1) 
and will be referred to as the open Friedmann separa¬ 
trix (OFS). The trajectories on the separatrix represent 
Milne models {x = 0, K = —1 and 0 ( 77 ) cx rj) which 


are equivalent to a Minkowski space-time in a hyper¬ 
bolic co-ordinate system. The second higher horizontal 
line {xA = — (a + 1)) is the separatrix which is the di¬ 
viding line between regions of phantom (x < xa) and 
non-phantom/standard behavior(a: > xa), we will call 
this the phantom separatrix (PS). The standard region 
corresponds to the case Bl, while the phantom region 
corresponds to the case Cl. In the phantom region the 
fluid violates the Null Energy Condition {p P < 0). 
This means the energy density is increasing (decreasing) 
in the future for an expanding (contracting) universe. In 
the standard case of the linear EoS in GR, this occurs 
when w < — 1 and ultimately leads to a Type I singu¬ 
larity [ 23 , E3' The parabola (y^ = x/3) represents the 
separatrix for flat Friedmann models {K = 0), we will 
call this the flat Friedmann separatrix (FFS). The inner 
most thick curve is the separatrix for closed Friedmann 
models {K = -|-1) and will be called the closed Friedmann 
separatrix (CFS). The separatrix has the form: 


y =3 


A(a l)x 
(a -I- 1 ) + tx_ 


3{o + l) 


(36) 


The constant A is fixed by ensuring that the saddle fixed 
point coincides with the fixed point representing the gen¬ 
eralized Einstein static model {E). The constant is given 
in terms of the EoS parameters and has the form: 


A= -- 


:(3q; -I- 1)(q; -|- 1) 


e(3a -I- 1) 
9 


3(0 + 1 ) 


(37) 


The Minkowski fixed point is located at the intersec¬ 
tion of the OFS and FFS. The generalized flat de Sit¬ 
ter fixed points are located at the intersection of the PS 
and FFS. The generalized Einstein static fixed point is 
located on the CFS. The trajectories between the OFS 
and the PS (0 < x < xa) represent models which exhibit 
phantom behavior (the case Cl). The open models in 
the phantom region are asymptotic to a Milne model in 
the past and to a generalized flat de Sitter model (dS'+) 
in the future. The closed models in the phantom re¬ 
gion evolve from a contracting de Sitter phase, through 
a phantom phase to an expanding de Sitter phase (phan¬ 
tom bounce). It is interesting to note that unlike the 
standard GR case the phantom behavior does not result 
in a Type I singularity but asymptotically evolves to a 
expanding de Sitter phase. This is similar to the behavior 
seen in the phantom generalized Chaplygin gas case . 
The trajectories on the PS all represent generalized de 
Sitter models (x' = 0). The fixed points represent gen¬ 
eralized flat de Sitter models [K = 0). The open model 
on the PS represent generalized open de Sitter models 
{K = — 1) in hyperbolic co-ordinates. The closed models 
on the PS evolve from a contracting phase to an expand¬ 
ing phase and represent generalized closed de Sitter mod¬ 
els {K = -|-1). The Friedmann equation can be solved for 
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FIG. 2: The normalized Hubble parameter, y for models with 
differing curvature. Starting from the top we have open, flat 
and the closed models. 



FIG. 3: The phase space for the high energy regime system 
with e = +1 and —1 < a < —1/3. The entire region corre¬ 
sponds to the case Al. 


such models to give: 

forA: = l, 

a{ri) = for A: = 0, (38) 

The region above the PS represents models which 
evolve in a non-phantom/standard manner (the case 
Bl). The trajectories in the expanding region {y > 0) of 
the phase space are asymptotic to a Type III singularity 
in the past®. The trajectories outside the FFS represent 
open models which evolve from a Type III singularity to 
a flat de Sitter phase, as do the trajectories on the FFS. 
The trajectories in between the CFS and the FFS evolve 
from a Type III singularity to a flat expanding de Sitter 
phase but may enter a phase of loitering. Loitering is 
characterized by the Hubble parameter dipping to a low 
value over a narrow red-shift range, followed by a rise 
again. In order to see this more clearly, we have plotted 
the normalized Hubble parameter (y) as a function of 
scale factor for three different trajectories in Fig. 2. The 
top two curves represent the open and flat models, with 


the Hubble parameter dropping off quicker for the flat 
Friedmann model. The lower most curve is the Hubble 
parameter for the closed model. The plot shows that the 
closed model evolves to a loitering phase. Loitering cos¬ 
mological models in standard cosmology were first found 
for closed FLRW models with a cosmological constant. 
More recently, brane-world models which loiter have been 
found , these models are spatially flat but can behave 
dynamically like a standard FLRW closed model. The in¬ 
teresting point here is that the models mentioned above 
loiter without the need of a cosmological constant (due to 
the appearance of an effective cosmological constant), the 
topology is asymptotically closed in the past and flat in 
the future. The trajectories inside the CFS can have two 
distinct types of behavior corresponding to the central 
regions above and below the generalized Einstein static 
fixed point. Trajectories in the lower region represent 
closed models which evolve from a contracting de Sitter 
phase, bounce and then evolve to a expanding de Sitter 
phase. The trajectories in the upper region evolve from 
a Type HI Singularity, expand to a maximum a (min¬ 
imum x) and then re-collapse to a Type HI singularity 
(we will refer to such re-collapsing models as turn-around 
models). 


The —l<a<—1/3 sub-case 

® The Type III singularity appears to be a generic feature of the 

high energy regime EoS and can occur both in the future and The phase Space for the System when —1 < O' < — 1/3 

the past. is shown in Fig. 3 . The fixed points representing the 
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FIG. 4: The phase space for the high energy regime system 
with e = +1 and a > —1/3. The entire region corresponds to 
the case Al. 


FIG. 5: The phase space for the high energy regime system 
with e = — 1 and a < — 1. The entire region corresponds to 
the case A2. 


flat generalized de Sitter models are no longer in the 
physical region {x > 0) of the phase space. The open, 
flat and closed Friedmann separatrices (OFS, FFS and 
CFS) remains the same. The phantom separatrix (PS) 
is no longer present and all trajectories represent mod¬ 
els with non-phantom/standard fluids (this corresponds 
to the case Al). The main difference is that the generic 
future attractor is now the Minkowski model. The trajec¬ 
tories between the OFS and FFS now evolve from a Type 
III singularity to a Minkowski model, as do the flat Fried¬ 
mann models. The models between the FFS and CFS 
now evolve from a Type III singularity to a Minkowski 
with the possibility of entering a loitering phase (as be¬ 
fore the model is asymptotically flat in the future). The 
trajectories inside the CFS and above the Einstein static 
fixed point still represent turn-around models. The tra¬ 
jectories inside the OFS and below the Einstein static 
model now represent standard bounce models, that is 
they evolve from a Minkowski model, contract to a finite 
size, bounce and then expand to a Minkowski model. 


3. The a > —1/3 sub-case 

Next we consider the system when a > —1/3, the 
phase space is shown in Fig. 4. The fixed point rep¬ 
resenting the Einstein static models is now located in the 
a; < 0 region of the phase space. The fluid in the entire 
physical region behaves in a non-phantom manner and 
corresponds to the case Al. The OFS and FFS remain 


the same and the CFS is no longer present. The trajec¬ 
tories between the OFS and FFS evolve from a Type III 
singularity to a Minkowski model. All trajectories above 
the FFS now represent turn-around models which start 
and terminate at a Type III singularity. The behavior of 
the models is qualitatively the same as that of the stan¬ 
dard FLRW model with a linear EoS where w = a, in 
the linear EoS case the Type III singularity is replaced 
by a standard “Big Bang”. 


D. The e = —1 case 


We now consider the system when we have a neg¬ 
ative quadratic energy density term (e = —1) in the 
high energy regime EoS. The character of the fixed point 
M is still undefined. The fixed point representing the 
generalized Einstein static model is now a center. The 
fixed points representing the expanding/contracting flat 
de Sitter points now have saddle stability. As before, 
the trajectories or fixed points in the y > 0 (y < 0) re¬ 
gion represent expanding models (contracting models), 
the black lines represent separatrix, grey lines represent 
example trajectories and fixed points are represented by 
black dots. 
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1. The a < —1 sub-case 


The phase space of the system when a < — 1 is shown 
in Fig. 5. The horizontal line {x = 0) is still the 
open Friedman separatrix (OFS). The parabola is the 
flat Friedmann separatrix (FFS). The intersection of the 
OFS and FFS coincides with the Minkowski fixed point. 
All the trajectories in the physical region of the phase 
space exhibit phantom behavior (corresponding to the 
case A2), the energy density increases in an expanding 
model. The trajectories in the expanding (contracting) 
region in general evolve to a Type III singularity in the 
future (past). The trajectories between the OFS and the 
FFS are asymptotic to a Milne model in the past and are 
asymptotic to a Type III singularity in the future. The 
trajectories on the FFS start from a Minkowski model 
and enter a phase of super-inflationary expansion and 
evolve to a Type III singularity. Trajectories that start 
in a contracting phase during which the energy density 
decreases, reach a minimum a (minimum x) and then ex¬ 
pand where the energy density increases represent phan¬ 
tom bounce models. The trajectories above the FFS rep¬ 
resent closed models which evolve through a phantom 
bounce, but start and terminate in a Type III singular¬ 
ity. The behavior of the models is qualitatively the same 
as that of the FLRW models with a phantom linear EoS 
{w < —1) except that the Type III singularity is replaced 
by a Type I (“Big Rip”) singularity. 





FIG. 6: The phase space for the high energy regime system 
with e = — 1 and — 1 < a < —1/3. The upper (lower) region 
corresponds to the case B2 (C2). 


2. The —l<a<—1/3 sub-case 

The phase space for the system when —1 < a < —1/3 
is shown in Fig. 6 . The lowest horizontal line (x = 0) 
is the OFS. The second higher horizontal line, xa = 
(a -T 1) is the phantom separatrix (PS), this divides the 
state space into regions of phantom (x > xa) and non¬ 
phantom/standard behavior (x < xa)- The phantom re¬ 
gion corresponds to the case B2 and the standard region 
corresponds to the case C2. The flat de Sitter (dS±) 
points are located at the intersection of the FFS and 
the PS. The open models in the standard matter region 
(0 < a: < Xa) are past asymptotic to open expanding de 
Sitter models in the past and evolve to Minkowski models 
in the future. The closed models in the region represent 
the standard bounce models, that is they evolve from a 
Minkowski model, contract to a minimum a (maximum 
x) and then expand to a Minkowski model. The trajec¬ 
tories above the PS (x > (a -I- 1)) all exhibit phantom 
behavior. The open models in this region are past asymp¬ 
totic to open de Sitter models in the past and evolve to a 
Type III singularity in the future. The closed models in 
the region all represent models which undergo a phantom 
bounce but start and terminate in a Type III singularity. 



1 


FIG. 7: The phase space for the high energy regime system 
with £ = — 1 and a > —1/3. The upper (lower) region corre¬ 
sponds to the case B2 (C2). 
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3. The a > —1/3 sub-case 

We now consider the system when a > —1/3, the phase 
space is shown in Fig. 7. The OFS, FFS and the PS are 
all still present, the phantom regions still corresponds to 
the case B2 and the standard region to the case C2. 
The trajectories in the phantom region {x > xa) behave 
in a similar manner to the previous case, as do the open 
models in the standard matter region (0 < x < xa). The 
main difference is in the region representing closed mod¬ 
els {K = 1) with non-phantom behavior. There is now 
a new fixed point which represents a generalized Ein¬ 
stein static model {E). The closed models in the region 
now represent oscillating models. This is represented by 
closed concentric loops centered on the Einstein static 
fixed point. These oscillating models also appear in the 
low energy system and will be discussed in more detail 
later. 


TABLE IV: Location and existence conditions of the fixed 
points of the low energy regime system. 

Name x y Existence 


(3q-I-1) 

0 

(30-1-1) 

ep 

1 / ~^p 

ep 

(a + l) 

' V 3(“+l) 

(o-t-l) 

ep 

/ -Ep 

ep 


V 3(c«+l) 

(o-l-l) 


point XA = Pk/\Po\ = 1), also xa > 0 for the 

points to be in the physical region of the phase space. 
The eigenvalues of the equilibrium points are given in 
Table El while the linear stability character is given in 

Table El 


IV. LOW ENERGY REGIME DYNAMICS 


A. The dimensionless dynamical system 


We now consider the system of equations for the low en¬ 
ergy regime EoS, which can be simplified and expressed 
in terms of the following dimensionless variables: 


The system of equations is then: 

2 X K 

y = 3"/W’ 

y' = ^ (3ep + (3a -h l)x), 

X = -3y (ep-I-(a-I-l)x). 


(39) 


(40) 

(41) 

(42) 


The discrete parameter Cp denotes the sign of the pres¬ 
sure term, Cp € { — 1,1}. The primes denote differentia¬ 
tion with respect to the new rj. The variables x and y are 
the new normalized energy density and Hubble parame¬ 
ter. As before only the positive energy density region of 
the phase space will be considered. The fixed points of 
the system and the existence conditions are given in Ta¬ 
ble HEI As before, by existence we mean the conditions 
on the parameters to insure x > 0 and x, y G R. 


The Minkowski model (x = y = 0) is no longer a fixed 
point of the system. The first fixed point (E) represents a 
generalized static Einstein model. This requires the over¬ 
all effective equation of state to be that of inflationary 
matter and therefore only exists when ep/(3a -I- 1) < 0. 
The last two points represent generalized expanding and 
contracting flat de Sitter models. These points only ex¬ 
ist if the fluid permits an effective cosmological constant 


TABLE V: Eigenvalues of the fixed points of the low energy 
regime system. 


Name 

Ai 

A 2 

E 



dS+ 

/-3(o-|-l) 

2 

V ^P 

-^-3ep(a + l) 

dS- 

/-3(0+1) 

2 

V -P 

-y/ —Sep (a + 1) 


TABLE VI: The linear stability of the fixed points for the 
low energy regime system. 


Name 

ep — +1 

ep — 1 

E 

Center {a yf —1/3) 

Saddle (a 7 ^ —1/3) 

dS+ 

Saddle (a < — 1) 

Attractor {a > —1) 

dS- 

Saddle (a < — 1) 

Repeller (a > — 1) 


B. The Cp = -1-1 case 

We start by considering the system when we have a 
positive constant pressure term (cp = -1-1) in the low en¬ 
ergy regime EoS. The Minkowski (x = y = 0) point is no 
longer present and the Einstein static solution has the 
stability character of a center. The fixed points repre¬ 
senting the generalized expanding/contracting de Sitter 
points {dS±) now have saddle stability. As before black 
lines represent separatrix, grey lines represent example 
trajectories and black dots represent fixed points of the 
system. 
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1. The a < —1 sub-case 



FIG. 8: The phase space for the low energy regime system 
with Ep = +1 and a < — 1. The upper (lower) region corre¬ 
sponds to the case El (FI). 


The phase space for the system when a < — 1 is shown 
in Fig. 8 . The open Friedmann separatrix (x = 0) 
is no longer present, and the x = y = 0 point is no 
longer a hxed point of the system. The horizontal line 
(xA = —(ct-l-l)”^) is the phantom separatrix (PS), divid¬ 
ing the state space into regions with phantom (x > xa) 
and standard {x < xa) behavior. The phantom region 
corresponds to the case El and the standard region cor¬ 
responds to the case FI. The parabola = x/3 is 
the separatrix representing the flat Friedmann models 
(FFS), this divides the remaining trajectories into open 
and closed models. The intersection of the PS and FFS 
coincides with the fixed points of the generalized flat de 
Sitter models. 

The trajectories in the upper region that start in a 
contracting phase (during which the energy density de¬ 
creases), reach a minimum a (minimum x) and then ex¬ 
pand, representing phantom bounce models which ter¬ 
minate in a Type I singularity. The closed models in 
the phantom region {x > x a) represent phantom bounce 
models which start and terminate in a Type I singularity 
The open models in the phantom region are asymp¬ 
totic to open de Sitter models in the past and evolve to 


a Type I singularity in the future. The trajectories be¬ 
low the PS {x < Xa) represent models which all behave 
in a standard manner (the FI case). The open models 
in this region are all non-physical as they all evolve to 
the X < 0 region of the phase space. The region cor¬ 
responding to closed models (above the FFS) contain 
a fixed point which represents the generalized Einstein 
static (E) model. The region is hlled by a infinite set 
of concentric closed loops centered on the Einstein static 
fixed point, the closed loops represent oscillating models. 
The Friedmann equation for such models is given by: 


2 / 


2 



Cp -h (a -h l)x 
B(a- 1 - 1 ) 


2 

3(ct + l) 


(43) 


The constant B is fixed by the location of the Einstein 
fixed point (B). The constant is given in terms of a and 


B = 


— 2 €„ 


(a -|- 1)(3 q; -|- 1) y — 


3o: -f 1 


3(0 + 1) 


(44) 


These oscillating models appear for —oo < a < —1/3 
when €p = -f-1 and are qualitatively similar to the oscil¬ 
lating models seen in the high energy case. The exact 
behavior of the variables for these models can be calcu¬ 
lated by hxing the EoS parameter a. The qualitative 
behavior remains the same for the models, however the 
maximum and minimum values of the variables change. 
In the case of a = —2/3 the equations are greatly sim¬ 
plified, the scale factor oscillates such that: 


i(r/) = tto (^1 + Vl - B sin(? 7 o - 77 )^ 


(45) 


The maximum and minimum scale factor is then: 


^max — {i+^/T^) 

j ^min — (i-VT^n^). (46) 

The normalized bubble parameter (y) is: 

\/l — K cos{r]o — v) 


y = Vo- 


1 -I- ^/l- K sin(77o - 77 ) 
The maximum and minimum y is given by: 


Vruax — Vo 


1 - K 


K 


ymin — yo 


l-K 


K 


The normalized energy density (x) is given by: 

1 — Vl — B sm(r/o — y)' 


, 1 -I- sin(r 7 o - r?) 

The maximum and minimum x are: 


(47) 


(48) 


(49) 


The Type I singularity is a generic feature of the low energy 
regime EoS and can appear both in the future and the past. 


fl + VT^\ fl-VT^\ 

( 50 ) 
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FIG. 9: The phase space for the low energy regime system 
with Sp = +1 and —1 < a < —1/3. The entire region corre¬ 
sponds to the case Dl. 


2. The —l<a<—1/3 sub-case 

We now consider the case when —1 < a < —1/3, the 
phase space is shown in Fig. 9. All trajectories in the 
physical region of the phase space exhibit standard be¬ 
havior and correspond to the case Dl. There is only one 
fixed point in the x > 0 region of the phase space, this 
point represents the generalized Einstein static model 
(E). The FFS represent models which evolve from a 
standard “Big Bang”, evolve to a Minkowski model and 
then to a standard “Big Crunch” (turn around model). 
The open models (below the FFS) are non-physical as 
the evolve into the x < 0 region. The trajectories above 
the separatrix represent closed models (K > 0) which 
oscillate indefinitely between a maximum and minimum 
a (minimum and maximum x), as seen in the previous 
case. 


3. The —1/3 < OL sub-case 

The phase space for the system when —1/3 < a is 
shown in Fig. 10. As in the previous sub-case the fluid 
behaves in a standard manner and corresponds to the 
case Dl. There are no fixed points in the physical re¬ 
gion of the phase space. The parabola is the FFS and 
represents flat model which evolve from a “Big Bang”, 
approach a Minkowski model and then re-collapse (turn¬ 
around models) to a “Big Crunch”. The open models 
(below the FFS) are all non-physical as they evolve to 



FIG. 10: The phase space for the low energy regime system 
with Cp = -1-1 and —1/3 < a. The entire region corresponds 
to the case Dl. 


the negative energy density region (x < 0) of the phase 
space. The closed models evolve from a “Big Bang”, 
reach a maximum a (minimum x) and re-collapse to a 
“Big Crunch”. 


C. The tp = —1 case 

We now consider the system when we have a nega¬ 
tive constant pressure term (cp = — 1) in the low energy 
regime EoS. As before,the Minkowski {x = y = 0) point 
is no longer a fixed point of the system and the OFS is 
not present. The fixed point representing the generalized 
Einstein static model (E) has the stability character of 
a saddle. The fixed points representing the generalized 
expanding (contracting) flat de Sitter points now have 
attractor (repeller) stability. 


1. The a < —1 sub-case 

The phase space for the low energy system when 
a < — 1 is shown in Fig. 11. All the trajectories in the 
X > 0 region of the phase space now exhibit phantom be¬ 
havior and correspond to the case D2. The open models 
are all non-physical as they all evolve from the negative 
energy density region of the phase space. The flat and 
closed models represent phantom bounce models which 
start and end in a Type I singularity. They evolve from 
a Type I singularity, contract, bounce at a minimum a 
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FIG. 11: The phase space for the low energy regime system 
with €p = —1 and a < —1. The entire region corresponds to 
the case D2. 



FIG. 12: The phase space for the low energy regime system 
with £p = —1 and — 1 < a < —1/3. The upper (lower) region 
corresponds to the case E2 (F2). 


(minimum x) and expand to a Type I singularity. 


S. The —1/3 < a sub-case 


2. The —l<a<—1/3 sub-case 

The phase space for the system when —1 < a < —1/3 
is shown in Fig. 12. The horizontal line, xa = (a + 1)”^ 
is the phantom separatrix (PS), dividing the state space 
into regions with phantom (x < xa) and standard be¬ 
havior (x > Xa)- The standard region corresponds to the 
case E2 and the phantom region corresponds to the case 
F2. The intersection of the PS and FFS coincides with 
the fixed points of the generalized flat de Sitter models 
(dS±). The flat expanding (contracting) de Sitter model 
is the generic future attractor (repeller). The open mod¬ 
els in the standard matter region (x > x a) evolve from a 
standard “Big Bang” to a flat expanding de Sitter phase. 
The closed models in this region evolve from a contract¬ 
ing flat de Sitter phase, reach a minimum a (maximum 
x), bounce and then evolve to expanding flat de Sitter 
phase. These models represent standard bounce models 
with asymptotic de Sitter behavior. The open models in 
the phantom region (x < xa) are all non-physical. The 
flat and closed models in this region represent models ex¬ 
hibiting phantom bounce behavior which avoid the “Big 
Rip” and instead evolve to an expanding flat de Sitter 
phase. 



FIG. 13: The phase space for the low energy regime system 
with Cp = — 1 and —1/3 < a. The upper (lower) region 
corresponds to the case E2 (F2). 
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We now consider the system in the parameter range 
— 1/3 < a, the phase space is shown in Fig. 13. The PS 
(Ja = (a + 1)“^), FFS = x/3) and generalized flat 
de Sitter points {dS±) still remain. The flat expanding 
(contracting) de Sitter model is the generic future attrac¬ 
tor (repeller). The inner most black curve is the closed 
Friedmann separatrix (CFS) and coincides with the gen¬ 
eralized Einstein static fixed point (E), which has saddle 
stability. The CFS is given by: 


y 


2 



{a -\- l)a: — 1 
2 


2 

3(a + l) 


(51) 


The constant H is a constant of integration and can be 
fixed by the location of the fixed point E. The constant 
is given in terms of a and has the form: 

Zl = (3a-k 1)”^^ . (52) 

The region below the PS {x < xa) remains qualita¬ 
tively the same. The open models in the standard mat¬ 
ter region (x > xa) all evolve from a “Big Bang” to a 
expanding flat de Sitter phase. The trajectories between 
the FFS and the CFS also evolve from a “Big Bang” to a 
generalized expanding flat de Sitter model with the pos¬ 
sibility of entering a loitering phase. The models inside 
the CFS can behave in one of two ways. The trajecto¬ 
ries above the generalized Einstein static point represent 
turn-around models which evolve from a “Big Bang”, 
reach a maximum a (minimum x) and then re-collapse 
to a “Big Crunch”. The trajectories below evolve from 
a contracting de Sitter phase to an expanding de Sitter 
phase and represent bounce models. 


V. THE FULL SYSTEM 


A. The dimensionless dynamical system 


We now consider the system of equations with the full 
quadratic EoS, this can be simplified in a similar fashion 
to the previous case by introducing a new set variables: 



\AfE\ ’ 


V = \/\fE\t , 


Po 

\/\pA 


(53) 


The system of equations then become: 

2 x K 
^ 3 \PcW 

y' = —y^ — - iZv + (3q; -I- l)a; -I- 3ex^), 
6 

x' = —iy {y -\- {a + l)a; -I- ex^) . 


(54) 

(55) 

(56) 


The parameter e denotes the sign of the quadratic term, 
e G { — 1, !}■ The parameter v is the normalized constant 
pressure term. The primes denote differentiation with re¬ 
spect to the new normalized time variable y and only the 


physical region of the phase space is considered {x > 0). 
The fixed points and their existence conditions are given 
in Table Eul The phase space undergo’s a topological 
change for special values of the ly parameter, these values 
can be expressed in terms of a and are: 


1^1 


(3a + 1)^ 
36 


1^2 


(a + 1)^ 
4 


(57) 


As in the previous case, by existence we mean a:: > 0 and 
x,y € M.. The general eigenvalues derived from the linear 
stability analysis are given in Table nsi The linear sta¬ 
bility character of the fixed points is given in Table 
The system has six fixed points and the sign of e no 
longer effects the linear stability character of the fixed 
point (but changes it’s position in the x — y plane). The 
first fixed point M represents a Minkowski model and 
is only present if = 0, the linear stability character is 
undefined. The second fixed point Ei represents an Ein¬ 
stein static model and has the linear stability character 
of a saddle. The third fixed point E 2 represents an Ein¬ 
stein static model with the linear stability character of a 
center. In general, this fixed point is surrounded by a set 
of closed concentric loops representing oscillating models. 
The next pair of fixed points dSi^± represents a set of gen¬ 
eralized flat de Sitter models, the expanding (contract¬ 
ing) model has attractor (repeller) stability. The next 
pair of fixed points dS 2 ,± also represents a set of general¬ 
ized flat de Sitter models, but now they have saddle sta¬ 
bility. The separatrix for open Friedmann models (OFS) 
is only present if = 0. The parabola y^ = x/3 (FFS) 
still separates the open and closed models. The separa¬ 
trix for the closed Friedmann models (CFS) is present 
for a narrow range of the parameters and always coin¬ 
cides with the fixed points representing the generalized 
Einstein static model, Ei. The fluid permits two possi¬ 
ble effective cosmological constants points, they are given 
by: 


PA,1 

(a + 1) 


IpcI 

2e 

+ 2e ’ 

Pa ,2 _ 

(a + 1) 


\Pc\ 

2e 

2e ■ 


(58) 

(59) 


Where <5 = (a -I- 1)^ — Aev. There is also a separatrix 
associated with each of the effective cosmological points, 
which divide the regions of phantom and non-phantom 
behavior. The separatrices will be referred to as the 
phantom separatrix (PSi which corresponds to the line 
X = XA,i)^ with the appropriate subscript. For special 
choices of parameters the separatrices coincide. The dis¬ 
cussion of the system will be split into the two categories, 
e = -1-1 and e = —1. 
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TABLE VII: The locations of the Hxed points of the fnll system. The existence conditions are also given, that is x,y,£ R 
and X > 0. To simplify the expressions special values of v are used which can be expressed in terms of a, these values are 




(3a+l)^ 


36 


and i /2 = 


Name 


X 


y 


Existence (e = +1) Existence (e = —1) 


M 


0 


0 




u = 0 


Si 


E 2 


dSi,± 


dS2,± 


(3a + l) 


6e 

' 6€ 

(3a + l) 

yj (Sa + l)^ —36eiv 

6e 

6€ 

(o+i) 

, •^(o' + l)2—4ez/ 

2e 

' 2e 

(c + l) 

yj (q: + 1)2 —4ei/' 

2e 

2e 


0 

0 


j- ( _ hi±i) V (a + l)'^-4ev 

* 6e 6e 


± - 


(g + l) _ V (a + l)'^-4eu 


u < ui, a < —1/3 
v < 0, a > —1/3 

0 < u < ui, a < —1/3 

u < U 2 , a < —1 
u < 0, a > —1 

0 < u < U 2 , a < —1 


—ui < u < 0, a > —1/3 


12 > 0 ,a<—1/3 
12 > — 121 , a > —1/3 


— V 2 < 12 < 0, a > —1 


12 > 0, a < —1 
—122 < t 2 , a > —1 


TABLE VIII: The linear stability character of the Hxed points 
for the full system. The stability character is only valid for 
choices of parameters which are consistent with the existence 
conditions and constraints given below. 


Name 

e = ±1 

Exceptions 

M 

Undefined 


El 

Saddle 

36ei2 7 ^ (3a -t 1)^ 

E 2 

Center 

36£U ^ (3a -h if 

dSi,+ 

Attractor 

4ei/ 7 ^ (a - 1 - 1 )^ 

dSi,- 

Repeller 

4ei2 7 ^ (a - 1 - If 

dS2,+ 

Saddle 

4ei2 7 ^ (q -I- if 

dS2,- 

Saddle 

4ei2 7 ^ (a -I- 1)^ 


B. The £ = +1 case 

We first consider the system when we have a positive 
quadratic energy density term (e = +1). The dynamical 
system can be further sub-divided into sub-cases with 
different values of the parameters a and v. The various 
subcases have been highlighted in Table HI 


The majority of sub-cases result in a phase space di¬ 
agram which is qualitatively similar to cases discussed 
in previous sections. That is the qualitative behavior of 
trajectories is the same even though the functional form 
of p{a) is different. The figure numbers not in bold (stan¬ 
dard text) indicate choices of variable for which the phase 
space is qualitatively similar to a previous case, with the 
following differences: 


• The regions which corresponded to different types 
of behavior of the fluid now change (replaced by 
new p{a) behavior): 

— The case D1 ^ Gl, 

— The case E2 ^ 13, 

— The case F2 —> 12, 

• The Type I singularities are now replaced by Type 
III singularities. 

There is a narrow range of the parameters for which 
the state space is qualitatively different. The figure num¬ 
bers given in bold in Tabled indicate the choice of vari¬ 
ables for which the state space is qualitatively different 
to previously discussed cases. We will now discuss the 
four sub-cases which are different to those discussed in 
previous sections. 


1. The a < —1, 12 — 121 sub-case 

The phase space of the system when a < — 1 and 
12 = i 2 i is shown in Fig. 14. As before the black lines 
represent separatrix, grey lines represent example trajec¬ 
tories and fixed points are represented by black dots. The 
fluid behaves in a standard manner and corresponds to 
the case Gl. This choice of parameters results in the two 
Einstein points (Ei) coinciding. The resulting fixed point 
is highly non-linear and cannot be classified into the stan¬ 
dard linear stability categories as in previous cases. The 
fixed point coincides with the CFS and the parabola is 
the FFS. The open models are all non-physical as they 
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TABLE IX: The Eigenvalues derived from the linear stability analysis of the fixed points for the full system. In order to simplify 
the form of the eigenvalues we introduce the following variables, 71 = (3a + 1), 72 = (a + 1) and 5 = ^(a + 1)^ — Aeiy. These 
eigenvalues are only valid for the choice of parameters consistent with the existence conditions. 

Name Ai A 2 


M 


El 


+ 


7 ^ — 71 'y'f —36ei/ — 36€i/ 

IsJ 




7^ —71 7^ — 36 €i^ —SGet' 

18^ 


E 2 


+ 


7^+71 Y/7i^-36eI7—36ei/^ ' 

18^ j 


7 j +71 ^ — Zdeiy—36ev' 

18^ 


dS2,± 


V^(i+¥) + ( 


6<5^(35-372-4)+8(72(3<5-1)+5) \ 2 
48€ 


V 


-(<^+ 72 ) 

6 e V 2 


-i)+ - 


66 ^( 36 + 372 + 4 )+ 8 ( 72 ( 33 + 1 )+ 3 ) \ 2 
48€ ) 


=y^(i + f)-( 


6<5^(35-372-4)+8(72(35-1)+5) \ 2 
48€ 




6i^(3a+372+4)+8(72(3^+l)+^) t 2 
48e / 


TABLE X: The various sub-cases of the full system when 
e = -fl. The figure numbers given in bold, indicate the choice 
of variables for which the state space is qualitatively different 
to previously discussed cases. 



a < —1 

-1 < a < -1/3 

-1/3 < a 

V > Vl 

FIG.IO 

FIG.IO 

FIG.IO 

U = Pi 

FIG.14 

FIG.14 

FIG.IO 

V2 < V < Vl 

FIG.15 

FIG.15 

FIG.IO 

V — V2 

FIG.16 

FIG.15 

FIG.IO 

0 < V < U2 

FIG.17 

FIG.15 

FIG.IO 

1 / = 0 

FIG.l 

FIG.3 

FIG.4 

u <0 

FIG.13 

FIG.13 

FIG.13 


evolve to the a; < 0 region of the phase space. The mod¬ 
els between the FFS and the CFS represent turn-around 
models which evolve from a Type III singularity^, evolve 
to a maximum a (minimum x) and then re-collapse. The 
trajectories above the CFS also represent similar turn¬ 
around models. 


2. The a < —1, V 2 < V < vi sub-case 

The phase space of the system when a < — 1 and 
r '2 < < i/i is shown in Fig. 15. The fluid behaves 

in a standard manner and corresponds to the case Gl. 
The Einstein hxed point of the previous case splits into 
two individual Einstein hxed points (Ei) via bifurcation. 
The hrst Einstein hxed point (Ei) coincides with the 


^ As with the case of the high energy EoS, the Type III singularity 
is a generic feature of the fully quadratic EoS. 



FIG. 14: The phase space for the full system with e = -|-1, 
a < — 1 and v — ui (additionally when a < —1/3 and u — vi). 
The entire region corresponds to the case Gl. 


CFS, while the second Einstein hxed point {E 2 ) is lo¬ 
cated inside the lower region enclosed by the CFS. Only 
the trajectories above the FFS differ from the previous 
case. The trajectories between the CFS and FFS still 
represent turn-around models which evolve from a Type 
III singularity but may now enter a loitering phase. The 
trajectories inside the CFS and above the Einstein hxed 
point (El), evolve from a Type III singularity, reach a 
maximum a and then re-collapse to a Type III singular- 
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FIG. 15: The phase space for the full system with e = +1, 
a < —1 and V 2 < v < (additionally when —l<a<—1/3 
and 0 < v < ui). The entire region corresponds to the case 

Gl. 


ity. The trajectories below Ei represent closed oscillating 
models, the closed loops are centered on the second Ein¬ 
stein fixed point (E 2 ). 


3. The a < — 1, v = V 2 sub-case 

The phase space of the system when a < —1 and 
V = V 2 IS shown in Fig. 16. The fluid behaves in a stan¬ 
dard manner in both regions and the upper (lower) region 
corresponds to the case H2 (HI). This choice of param¬ 
eters results in the two sets of generalized de Sitter points 
(dSi^±) coinciding. The fixed points coincide with PSi 
{x = XA^i) and the FFS. The resulting fixed points are 
highly non-linear, the points have shunt stability along 
the FFS direction and the generalized expanding (con¬ 
tracting) de Sitter point has attractor (repeller) stability 
along the PSi direction. The two Einstein points (Ei) 
and the CFS are still present. In the x < XA,i region, 
the open models are all non-physical as they evolve to 
the X < 0 region and the closed models represent oscil¬ 
lating models which are centered on the Einstein point 
{E 2 ) with center linear stability. In the x > XA,i region, 
the open models are asymptotic to a Type III singularity 
in the past and a expanding flat de Sitter phase {dSi^+) 
in the future. The trajectories between the FFS and the 
CFS evolve from a Type III singularity to dSi^+ with 
the possibility of entering a loitering phase. The models 
inside the CFS and above the Ei point represent turn¬ 



-1 -0.5 0 0.5 1 
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FIG. 16: The phase space for the full system with e = -|-1, 
a < — 1 and v = V 2 - The upper (lower) region corresponds to 
the case H2 (HI). 


around models which asymptotically approach a Type III 
singularity. The closed models inside the CFS and below 
the E 2 point are asymptotic to a contracting de Sitter 
model phase (dSi^-) in the past and a expanding de Sit¬ 
ter phase (dSi^+) in the future. The generic attractor in 
the X > XA,i region is the dSi^+ fixed point. 


4- The a < —1, 0 < v < U 2 sub-case 

The phase space of the system when a < — 1 and 
0 < ^ < r '2 is shown in Fig. 17. The upper (lower) 
horizontal line is the (F’S' 2 ). The region above PSi 
corresponds to the case 13 and is qualitatively similar to 
the H2 region in the previous sub-case. The region be¬ 
low PS 2 corresponds to the case II and is qualitatively 
similar to the HI region in the previous sub-case. The 
set of generalized fiat de Sitter fixed points {dSi^±) of 
the previous case split into two sets of generalized flat 
de Sitter fixed points via bifurcation. The upper (lower) 
set of generalized de Sitter points, dSi^± {dS 2 ,±) have 
attractor/repeller (saddle) stability. The region between 
PSi and PS 2 corresponds to the case 12 and the fluid 
behaves in a phantom manner. The open models in this 
region are asymptotic to open de Sitter models in the 
past and flat de Sitter models in the future. The closed 
models in the phantom region represent phantom bounce 
models which asymptotically approach a expanding (con¬ 
tracting) de Sitter phases in the future (past). 
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FIG. 17: The phase space for the full system with e = +1, 
a < —1 and Q < v < V 2 - The upper, middle and lower regions 
correspond to the cases 13, 12 and II respectively. 


FIG. 18: The phase space for the full system with e = —1, 
a > —1/3 and —ui < v < G. The upper, middle and lower 
regions correspond to the cases 16, 15 and 14 respectively. 


C. The e = —1 case 

We now consider the system when we have a negative 
quadratic energy density term (e = — 1). As before the 
system can be sub-divided into various sub-cases with 
different values of parameters of a and v. The various 
sub-cases have been highlighted in Ta.ble fxTI 


TABLE XI: The various sub-cases of the e = — 1 full system. 
The hgure numbers given in bold, indicate the choice of vari¬ 
ables for which the phase space is qualitatively different to 
previous cases. 



a < —1 

-1 < a < -1/3 

-1/3 < a 

ly >0 

FIG.8 

FIG.8 

FIG.8 

iy = 0 

FIG.5 

FIG.6 

FIG.7 

—i/1 < i/ < 0 

FIG.ll 

FIG.20 

FIG.18 

ly = —i/i 

FIG.ll 

FIG.20 

FIG.19 

— i/2 < i/ < —i/1 

FIG.ll 

FIG.20 

FIG.20 

ly = —V2 

FIG.ll 

FIG.21 

FIG.21 

V < —i/2 

FIG.ll 

FIG.ll 

FIG.ll 


As before, the figure numbers not in bold (standard text) 
indicate choices of variable for which the phase space is 
qualitatively similar to a previous case, with the following 
differences: 


• The regions which corresponded to different types 
of behavior of the fluid now change (replaced by 
new form of pifl))'. 

— The case D2 ^ G2, 

— The case El 16, 

— The case FI ^ 15, 

• The Type I singularities are now replaced by Type 
III singularities. 

There are choices of parameters for which the phase space 
is different (figure numbers in bold in Tabie lxTll and these 
four sub-cases will be discussed in the following sections. 


1. The a > —1/3, —v\ < v < 0 sub-case 

The phase space of the system when a > —1/3 and 
—vi < < 0 is shown in Fig. 18. The upper (lower) 

horizontal line x = {x = xa.i) is the PS2 (PSi) 
(they have swapped position with respect to the e = -1-1 
case). The region above PS 2 corresponds to the case 16 , 
the region below PSi corresponds to the case 14 and the 
fluid behaves in a phantom manner in both regions. The 
region between PSi and PS 2 corresponds to the case 15 
and the fluid behaves in a standard manner. The lower 
set of generalized de Sitter points {dSi^± - at the inter¬ 
section of PSi and FFS) have attractor/repeller stability, 
while the upper set {dS 2 ,± - at the intersection of PS 2 
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FIG. 19: The phase space for the full system with e = — 1, 
a > —1/3 and v = —vi. The upper, middle and lower regions 
correspond to the cases 16, 15 and 14 respectively. 


and FFS) have saddle stability. The CFS is located in be¬ 
tween PSi and PS 2 and coincides with the Einstein point 
{El). The open models in the x < xtii region (the case 
14) are all non-physical as they evolve from the a: < 0 re¬ 
gion of the phase space. The closed models in this region 
represent phantom bounce models which evolve from a 
contracting de Sitter phase {dSi^-) to a expanding de 
Sitter phase {dSi^+). The open models in the standard 
region (xAp < x < xa ,2 corresponding to the case 15 ) are 
asymptotic to a generalized open de Sitter model in the 
past and generalized flat de Sitter model in the future 
(the future attractor has lower x and y). The models 
between the CFS and the FFS in this region represent 
bounce models which evolve from a contracting de Sitter 
phase to a expanding de Sitter phase with the possibility 
of entering a loitering phase. The models enclosed by the 
CFS can be split into two groups. The models above the 
fixed point, Ei represent oscillating models, the closed 
loops are centered on the fixed point E 2 . The models 
below the fixed point, Ei represent bounce models which 
evolve from dSi,- to dSi^+. In the x > xa,2 region (the 
case 16) the open models are asymptotic to generalized 
open de Sitter models in the past and a Type III sin¬ 
gularity in the future. The closed models in this region 
represent phantom bounce models which evolve from a 
Type III singularity, reach a minimum a (minimum x) 
and then evolve to a Type III singularity. The gener¬ 
alized expanding flat de Sitter model, dSi^+ (Type III 
singularity) is the generic future attractor in the region 
X < XA .2 {x > XA. 2 )- The trajectories in the regions, 



FIG. 20: The phase space for the full system with e = — 1, 
a > —1/3 and —V 2 < v < —vi (additionally when —1 < a < 
— 1/3 and — V 2 <v<(^). The upper, middle and lower regions 
correspond to the cases 16, 15 and 14 respectively. 


X < XA,i and x > xa,2 remain qualitatively similar in the 
following two cases (Fig.19, 20). 


2. The a > —1/3, v = —vi sub-case 

The phase space of the system when a > —1/3 and 
v = —vi is shown in Fig. 19. The phase space is 
equivalent to the previous sub-case, except for the re¬ 
gion xa,i < X < xa ,2 (the case 15). The open models 
in this region are still asymptotic to generalized open 
(flat) de Sitter models in the past (future). The behav¬ 
ior of the closed models has now changed, there are no 
longer trajectories representing oscillating models. The 
two generalized Einstein fixed points (Ei) have now coa¬ 
lesced to form one fixed point via bifurcation. The closed 
models above Ei represent bounce models which evolve 
from dSi^- to dSi,+ , with the possibility of entering a 
loitering phase. The closed models below Ei represent 
bounce models which evolve from dSi^- to with¬ 

out entering a loitering phase. 


3. The a > —1/3, —122 < v < —ui sub-case 

The phase space of the system when a > —1/3 and 
— 1^2 < V < —vi is shown in Fig. 20. The phase space 
is qualitatively similar to the previous sub-cases except 
for the XA,i < x < xa,2 region. There are no longer 
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VI. DISCUSSION AND CONCLUSIONS 



FIG. 21: The phase space for the full system with e = — 1, 
a > —1/3 and v = —V 2 (additionally when —1 < a < —1/3 
and V = —V 2 )- The upper (lower) region corresponds to the 
case H4 (H3). 


any fixed points representing generalized Einstein static 
models and the CFS is no longer present. The open 
models in the region behave as in previous sub-cases. 
The closed models in the region represent bounce model, 
which evolve to a expanding (collapsing) de Sitter phase 
in the future (past) without the possibility of entering a 
loitering phase. 


/. The a > —1/3, v = —V 2 sub-case 

The next case is the phase space of the system when 
a > —1/3 and ly = —V 2 and is shown in Fig. 21. The fluid 
behaves in a phantom manner in both regions and the up¬ 
per (lower) region corresponds to the case H4 (H3). The 
two sets of generalized de Sitter points (dS'i,±) have now 
coalesced into a single set of generalized de Sitter points 
{dS±) which are located at the intersection of the FFS 
and the PSi which have also coalesced to form a single 
separatrix (xa^ = The resulting fixed points are 

highly non-linear, the points have shunt stability along 
the FFS direction and the generalized expanding (con¬ 
tracting) de Sitter point has attractor (repeller) stabil¬ 
ity along the PSi direction. The Type III singularity 
is the generic attractor in the upper region {x > x\^i) 
and the dSi^+ is the generic attractor in the lower region 
{x < XAy). 


In this paper we have systematically studied the dy¬ 
namics of homogeneous and isotropic cosmological mod¬ 
els containing a fluid with a quadratic EoS. This has it’s 
own specific interest (see Section I for a variety of moti¬ 
vations) and serves as a simple example of more general 
EoS’s. It can also be taken to represent the truncated 
Taylor expansion of any barotropic EoS, and as such it 
serves (with the right choice of parameters) as a useful 
phenomenological model for dark energy, or even UDM. 
Indeed, we have shown the dynamics to be very different 
and much richer than the standard linear EoS case, find¬ 
ing that an almost generic feature of the evolution is the 
existence of an accelerated phase, most often asymptot¬ 
ically de Sitter, thanks to the appearance of an effective 
cosmological constant. Of course to properly build phys¬ 
ical cosmological models would require to consider the 
quadratic EoS for dark energy or UDM together with 
standard matter and radiation. Our analysis was aimed 
instead to derive and classify the large variety of different 
dynamical effects that the quadratic EoS fluid has when 
is the dominant component. In this respect, it should 
be noticed that a positive quadratic term in the EoS al¬ 
lows, in presence of another fluid such as radiation, equi- 
density between the two fluid to occur twice, i.e. the 
quadratic EoS fluid can be dominant at early and late 
times, and subdominant in an intermediate era. 

In Section II we have made some general remarks, 
mostly based on conservation of energy only and as such 
valid independently of any specific theory of gravity. We 
have also given the various possible functional forms of 
the energy density as a function of the scale factor, p{a ), 
and listed the many subcases, grouped in three main 
cases, what we call: i) the high energy models (no con¬ 
stant Po term); ii) the low energy affine EoS with no 
quadratic term; Hi) the complete quadratic EoS. 

The quadratic term in the EoS affects the high energy 
behavior as expected but can additionally affect the dy¬ 
namics at relatively low energies. First, in Section III, we 
have concentrated on the high energy models. The spe¬ 
cific choice of parameters fixes the behavior of the fluid, 
it can behave in a phantom or standard manner. In the 
case of phantom behavior, p can tend to zero at early 
times and either tend to an effective cosmological con¬ 
stant (Cl) or a Type III singularity (A2) at late times. 
Alternatively p can also tend to an effective cosmological 
constant in the past (B2) and a Type III singularity at 
late times. When the fluid behaves in a standard man¬ 
ner, it can tend to a Type III singularity at early times, 
with p either tending to zero (Al) or to an effective cos¬ 
mological constant (Bl) at late times. The fluid can 
also behave as an effective cosmological constant at early 
times with p decaying away at late times (C2). The ef¬ 
fective cosmological constant allows for the existence of 
generalized Einstein static {E) and flat de Sitter fixed 
{dS±) points which modify the late time behavior. The 
main new feature is the existence of models which evolve 
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from a Type III singularity and asymptotically approach 
a flat de Sitter model (d5'+). Of specific interest are the 
closed models of this type, which can also evolve through 
an intermediate loitering phase. 

Neglecting the quadratic term, in Section IV we have 
considered the low energy models with affine EoS. As ex¬ 
pected, the constant term in the quadratic EoS affects 
the relatively low energy behavior. It can result in a va¬ 
riety of qualitatively different dynamics with respect to 
those of the linear EoS case. Again, the fluid can have a 
phantom or standard behavior. When the fluid behaves 
in a phantom manner, p can tend to an effective cosmo¬ 
logical constant (F2), or can tend to a Type I (“Big Rip”) 
singularity (D2) at late times. Alternatively, p can also 
tend to an effective cosmological constant in the past and 
a Big Rip in the future(El). When the fluid behaves in 
a standard manner, we recover the linear EoS at early 
times and p can either tend to zero (Dl) or to an effec¬ 
tive cosmological constant (E2) at late times. The fluid 
can also behave as an effective cosmological constant at 
early times, with p decaying away at late times (FI). 
The effective cosmological constant allows for the exis¬ 
tence of new fixed points(if and dS±). Comparing with 
standard linear EoS cosmology, the most interesting dif¬ 
ferences are new closed models which oscillate indefinitely 
and new closed models which exhibit phantom behavior 
which do not terminate in a “Big Rip”, but asymptot¬ 
ically approach an expanding flat de Sitter model (flat 
and closed models where the fluid behaves as case F2). 

When we study the dynamics of the system with the 
complete quadratic EoS, Section V, we see the appear¬ 
ance of new fixed points representing generalized Ein¬ 
stein and de Sitter models which are not present in the 
high/low energy systems. The various models of the sim¬ 
plified systems are present in the full system (but with 
differing p(a)), but there are also models with qualita¬ 
tively new behavior. As with the previous cases, in the 
case of phantom behavior, p can tend to zero at early 
times and either tend to an effective cosmological con¬ 
stant (H3 and 14) or a Type III singularity (G2) at late 
times. Alternatively p can also tend to an effective cos¬ 
mological constant in the past (H4 and 16) and a Type 
III singularity at late times. Einally, in the phantom 
case p can also tend to an effective cosmological constant 
both in the past and future (12). In the case of standard 
behavior the fluid can tend to a Type III singularity at 
early times, with p either tending to zero (Gl) or to 


an effective cosmological constant (H2 and 13) at late 
times. The fluid can also behave as an effective cosmo¬ 
logical constant at early times with p decaying away at 
late times (HI and II). Finally, in the standard fluid 
case p can also tend to an effective cosmological constant 
both in the past and future (15). There are models which 
evolve from a Type III singularity, reach a maximum a 
(minimum x) and then evolve to Type III singularity. 
These also enter a loitering phase before and after the 
turn around point. We also see bounce models which 
enter a loitering phase and asymptotically tend to gen¬ 
eralized expanding (contracting) de Sitter models at late 
(early) times. 

Of specific interest are models which evolve from a 
Type III singularity as opposed to the standard “Big 
Bang” (Al, Bl). The simplest models of this type corre¬ 
spond to the high energy EoS with a positive quadratic 
term (is possible to recover standard behavior at late 
times). For these models the positive quadratic energy 
density term has the potential to force the initial sin¬ 
gularity to be isotropic. The effects of such a fluid on 
anisotropic Bianchi I and V models is investigated in Pa¬ 
per II |48l| . This is achieved by carrying out a dynamical 
systems analysis of these models. Additionally, using a 
linearized perturbative treatment we study the behav¬ 
ior of inhomogeneous and anisotropic perturbations at 
the singularity. The singularity is itself represented by 
an isotropic model and. If the perturbations of the latter 
decay in the past, this model represents the local past at¬ 
tractor in the larger phase space of inhomogeneous and 
anisotropic models (within the validity of the perturba¬ 
tive treatment). This would mean that in inhomogeneous 
anisotropic models with a positive non-linear term (at 
least quadratic) in the EoS isotropy is a natural outcome 
of generic initial conditions, unlike in the standard linear 
EoS case where generic cosmological models are, in GR, 
highly anisotropic in the past. 
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